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Abstract. We consider an extension of the lepto-specific 2HDM with an extra singlet S as 
a dark matter candidate. Taking into account theoretical and experimental constraints, we 
investigate the possibility to address both the 7 -ray excess detected at the Galactic centre 
and the discrepancy between the Standard Model prediction and experimental results of 
the anomalous magnetic moment of the muon. Our analyses reveal that the SS —> t + t~ 
and SS —> bb channels reproduce the Galactic centre excess, with an emerging dark matter 
candidate which complies with the bounds from direct detection experiments, measurements 
of the Higgs boson invisible decay width and observations of the dark matter relic abundance. 
Addressing the anomalous magnetic moment of the muon imposes further strong constraints 
on the model. Remarkably, under these conditions, the SS —> bb channel still allows for the 
fitting of the Galactic centre. We also comment on a scenario allowed by the model where the 
SS —> t + t~ and SS —> bb channels have comparable branching ratios, which possibly yield 
an improved fitting of the Galactic centre excess. 
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1 Introduction 

The matter content of the Universe is dominated by a weakly interacting component, the 
Dark Matter (DM), which accounts for about 26 % of the energy density of the former at 
the present time. According to the current paradigm, DM consists of a relic abundance of 
weakly-interacting massive particles (WIMPs). This scheme is motivated by the fact that 
particles with masses and annihilation cross sections set by the electroweak scale, a known 
scale of Nature, yield through the freeze-out mechanism 1 a relic abundance comparable with 
the observed DM one. Given that the freeze-out mechanism is a natural consequence of an ex¬ 
panding Universe and that new weakly-interacting particles arise in well motivated extensions 
of the Standard Model (SM), e.g. super-symmetric theories, the above observation is referred 
to in literature as the “WIMP miracle”. The WIMP paradigm motivates the experimental 
efforts behind different DM searches: for instance, the direct detection experiments rely on 
the possible elastic scattering between DM and Standard Model (SM) particles mediated by 
the weak force. Likewise, the weak interactions could allow for the DM production at collid¬ 
ers, motivating the dedicated collider searches. Finally, the annihilation of DM particles in 
space, mediated by the same force, provide the basis for DM indirect detection. Regarding 
the latter possibility, notice that DM annihilations take place only in regions where the DM 
abundance is sizeable and, therefore, possible annihilation signals are expected from dense 
DM regions as the Galactic Centre (GC) and the dwarf satellite galaxies of the Galaxy. In 
addition to that, DM annihilations in space yield large scale effects that can be investigated, 
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for example, by their impact on the Cosmic Microwave Background radiation. The latter is 
especially constricting for light WIMPs, < 10 GeV [3-10]. 

In this light, the DM halo at GC should then provide the strongest indirect detection 
signals and, interestingly, in 2009 a spatially extended 7 -ray excess was detected in this region 
by the Fermi LAT telescope in the energy range of 1-5 GeV [11]. We refer to this signal, con¬ 
firmed by many other studies [12-23], as the Galactic Centre Excess (GCE). Clearly the GC 
harbours an extremely rich environment populated with stars, stellar relics, cosmic rays, dust, 
gas and the central black hole and, so far, it has not been possible to clearly disentangle the 
potential annihilation signal from the partially known astrophysical background due, for in¬ 
stance, to millisecond pulsars [24] and ultra-energetic events form the past [25]. Nevertheless, 
the DM annihilation interpretation is supported by the estimates of the required annihilation 
cross section, matching the order of the thermal freeze-out one, as well as by the morphology 
of the signal region. 

When interpreting the detected GCE in the light of the WIMP paradigm it is necessary 
to address the strong constraints from the direct detection experiments XENONIOO [26] and 
LUX [27], which disfavour the characteristic weak-scale values of scattering cross section. This 
difficulty has pushed the scientific community to consider alternative models that possess an 
annihilation cross section large enough to explain the detected signal, but present a suppressed 
DM-nucleon scattering cross section. An example is provided by the so-called Coy Dark 
Matter [28] models, in which the DM particle is a Dirac fermion x interacting with the SM 
particles through a light (~10 GeV) pseudoscalar mediator. By assuming that the couplings to 
the SM particles are proportional to the corresponding Higgs Yukawa couplings, as motivated 
by minimal flavour violation [29], the model avoids the tight direct detection constraints 2 . 
Furthermore, if the DM particles are lighter than the top quark, the dominant annihilation 
channel is XX bb and with the choice m x 6 [40, 55] GeV, the GCE can be fitted for a 
natural value of the DM annihilation cross section [18-23]. 

Alternatively, it is possible to fit the signal with the XX —> rr channel by adopting 
a DM mass m x ~ 10 GeV. In this case, assuming leptophilic SM-pseudoscalar couplings 
proportional to the corresponding lepton Yukawa couplings also removes the tension between 
the SM prediction and the measurement of anomalous magnetic moment of the muon a^ [31]. 
In this regard, the long-standing ~ 3cr discrepancy 3 that afflicts this quantity can be regarded, 
along with the origins of DM, neutrino masses and baryon asymmetry of the Universe, as a 
compelling experimental evidence that points to physics beyond the SM. In this paper we 
adopt the above point of view and fit the GCE within the framework of the extended Two 
Higgs Doublet Model (2HDM), investigating also the consequences of the implied particle 
content on the anomalous magnetic moment of the muon. 

The 2HDM [34, 35] is a well known extension of the SM in which a second scalar 
doublet field is added to the particle content of the theory. The latter transforms under the 
symmetries of the SM in the same way as the original Higgs doublet and, depending on the 
model, is allowed to develop a non-zero vacuum expectation value (VEV). The SM-like Higgs 
boson, discovered at the LHC [36, 37], is a combination of the two doublets. Barring models 
that induce flavour-changing neutral currents either at the tree or loop level, there are four 
possible ways to assign the Yukawa couplings with the SM fermions to the two scalar doublets 
[35, 38, 39]. As shown in table 1, these correspond to different Z 2 parity assignments, which 

2 Interestingly, a light mediator can be exploited for the annual modulation signal of DAMA [30]. 

,s For a review see [32] and references therein. We consider a conservative value for the statistical significance 
of the mentioned discrepancy [33] . 
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mark the four possible 2HDM variants: type-I, II, X (or “lepto-specific”) and Y (or “flipped”). 
The model in which the additional Higgs doublet does not acquire a VEV nor couples to the 
fermions is instead called the inert doublet model [40, 41]. 

Table 1: Z 2 parity and Yukawa enhancement factors y l x . 
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In this paper we focus on an extension of the lepto-specific 2HDM, where a further real 
singlet scalar field S is added to play the role of dark matter [42-50] (see [51] for a recent 
review). Whereas such a construction has been previously considered in literature [52-54] 
(or as an effective field theory [55]), the novelty of our work is in the detailed analysis of the 
GCE fit, which accounts for the latest bounds that the mentioned direct detection experiments 
yield. We furthermore scrutinise the impact of the extended particle content presented by the 
model on the anomalous magnetic moments of the electron and muon, investigating whether 
the same values of the parameters that fit the GCE reduce, as well, the tension that afflict 
the SM prediction of the latter. 

The paper is organised as follows: in the next section we provide the details of the 
considered model, in order to address in section 3 the bounds that perturbativity, perturbative 
unitarity, vacuum stability and Electroweak precision data impose. Section 4 is dedicated to 
fitting the GCE, while in section 5 we discuss the muon anomalous magnetic moment in the 
adopted framework. We summarise our results in section 6 . 


2 The considered model 


The Higgs sector of the 2HDM presents two 51/(2) scalar doublets H\ and H 2 . In order 
to avoid flavour-changing neutral currents, it is customary to impose a Z 2 symmetry on the 
scalar potential. As said before, there are four possible ways to assign a Z 2 parity to the 
fields, summarised in table 1 . 

In the lepto-specific, or type-X 2HDM, H\ = Hi couples only to leptons while H 2 = Hq 
only to quarks. Our scalar sector contains, in addition to the two doublets, a real singlet S 
that is our dark matter candidate . 4 We forbid explicit CP-violation but allow for the soft 
breaking of the Z 2 symmetry that characterises the considered model. The most general 
scalar potential in line with these assumptions is then 


1 / = - ^ Q H Q - \ Q {H\H Q + H' Q H L ) + -y s S* + A S 5 4 

+ A 4H[H L ) 2 + A 2{H ] Q H Q f + A 3 (h[h l )(H^ q H q ) + A A {H[H Q ){H ] Q H L ) 
(- H' Q H L ) 2 + (h[Hq) 2 ] + A slH[h l S 2 + \sqH' q H q S 2 , 


( 2 . 1 ) 


+ ^A 5 


4 Similar models have been explored in the context of DM stabilised by Zjv symmetries in [56, 57]. 
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where both /jl 2 L q and A 5 are real and the masses // 2 q are responsible for the soft breaking of 
the Z 2 symmetry. 

As usual we parametrise the degrees of freedom contained in the Higgs doublets as 5 


H = ( h ‘ 

\(vi + hi — iai)/y/2 

where v\ and V 2 are the VEVs of the two scalar doublets. The particle content of the SM is 
then effectively extended to accommodate an extra neutral scalar H, the charged scalars H± 
and a neutral pseudoscalar A. We indicate with h the usual SM Higgs boson. 

In the following, we choose as free parameters the SM-like Higgs mass rrih = 125.09 GeV 
[59], the pseudoscalar mass m A , the neutral scalar mass mu, the charged Higgs mass mu±, 
the DM mass ms, Ai, A 2 , A = A 3 + A 4 +A 5 , A 5 , A sl and A sq, the neutral scalar mixing angle 
a, together with v 2 = v\ + v 2 = 246.2 GeV and tp = tan /3 = 1 ) 2 /v\. 

We can express 
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( 2 . 6 ) 
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where M 2 = p? L Q (tp + ^-j, sp- a = sin {/3 — a) and cp- a = cos {/3 — a) [58]. Because the 

coupling Ai, as given by eq. (2.3), tends to large for increasing tp, we adopt Ai as an input 
parameter instead of Plq, in virtue of eq. (2.3). 


3 Theoretical and experimental constraints 
3.1 Perturbativity and perturbative unitarity 

To satisfy a perturbativity constraint, we require the absolute values of all quartic couplings 
to be smaller than 47 T. 

We also consider perturbative unitarity. Because at high energy the scattering rates are 
dominated by quartic contact interactions, we require the eigenvalues of the s-wave 5-matrix 
of two-to-two scalar boson states to satisfy 

1 , X 

I a ij I V —. (3.1) 


5 We follow the conventions of [58]. 
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For the 2HDM the unitarity constraints have been given in [60-63]. As in [56], we extend 
the method of [62, 63] to our field content. The eigenvalues are given by 


167TAS2 = Ai + A 2 ± \J (Ai — A 2 ) 2 + A§, (3.2) 

16vrA^ d = A 3 + A 4 , (3.3) 

167rA^ d = A 3 - A 4 , (3.4) 

IOttAq^jI 1 = Ai + A 2 ± \J (Ai — A 2) 2 + A|, (3-5) 

16-7rAoo± = A 3 + 2A 4 ± 2A 5 , (3.6) 


and the Aqq^ 3 are given by x the three roots of the polynomial equation 

2.x 3 — x 2 (12Ai + I 2 A 2 + y/2\s) + I 2 A 2 A 52 , + (2A 3 + A 4 )(\/2(2A 3 + A 4 )As — 4Ag^A5Q) 

+ 12Ai(— 3 A/ 2 A 2 A 5 + A|q) — 2x((2A 3 + A 4 ) 2 — 3 ( 12 A 4 A 2 + y/2X\X g + \^2X2Xg) (3-7) 

+ tfsL + ^Sq) = 0 . 

3.2 Vacuum stability 

The potential is bounded below in the limit of large field values if the quartic couplings satisfy 
the copositivity conditions [64] 

Ai >0, A 2 > 0, A s > 0, (3-8) 

2\/AiA2 + A 3 + A 4 — |A 5 1 > 0 2\JX\Xs + Xsl > 0 2 \J A 2 A s + A sq > 0 (3.9) 

2 \J A 2 A 2 A s + V^A^q + \fX 2 XsL + \/As(A 3 + A 4 — | A 5 I) (3.10) 

+ \j (2\/AiA ~s + Xsl){^V^Xs + Xsq){2\/ A 1 A 2 + A 3 + A 4 — |As|) > 0. 

We have derived the renormalisation group equations (RGEs) for our model with the 
help of the PyR@TE package [65]. The /3-functions are given in the Appendix A. We use 
the RGEs to find the scale A at which the model becomes non-perturbative or loses vacuum 
stability. The results are shown in figure 1 as a function of tua and run- 


3.3 Electroweak precision data 

The electroweak oblique parameters S, T and U parametrise the effect of new physics on the 
gauge bosons two point functions [ 66 ]. 

The latest results from the GFitter group [67] are 

AS = 0.05 ± 0.11, AT = 0.09 ±0.13. (3.11) 


In the SM-like limit sin(/3 — a) —> 1, the electroweak precision S, T and U parameters 
[ 66 ] reduce to those of the inert doublet model [40]. If <C mz m#± ~ mu, they 
are [58] 


AS « 0 . 022 , 


\m H ± — m H 

AT « mH——^ - 2 “ 

327r z a em v^ 


(3.12) 


Since the T parameter constrains |mn — i^h ± \ = 0(10) GeV, in the following we always 
consider m#± = rriH to comply with this constraint. 
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Figure 1: The scale A at which the model becomes non-perturbative as a function of m#± = 
mu vs. m A . 

3.4 \ hAA coupling 

The SM-like Higgs boson can decay into two pseudoscalars A if m A < m^J 2. Dedicated 
collider experiments bound the corresponding branching ratio to be smaller than 0.2 [ 68 ], 
therefore we opt to set the coupling A hAA to zero [58] by imposing that 



(3.13) 


With the above prescription, the model still recovers the SM-like behaviour of h in 
the limit tp 1, allowing however for testable consequences on sign of the lepton Yukawa 
couplings [58]. 

4 Fitting the 7 -ray excess at the Galactic Centre 
4.1 The 7 -ray excess at the Galactic Centre 

The data provided by the Fermi LAT 7 -ray telescope [11] show a spatially extended 7 -ray 
excess in the range 1-5 GeV at the GC [12-23]. In this paper we assume that the signal 
originates in DM annihilation, neglecting possible competing astrophysical processes as mil¬ 
lisecond pulsars [24] and past energetic events in the GC [25]. If the GCE is fitted by using 
a single channel, i.e. a single SM final state, the best fit is provided by bb and cc. Lighter 
or heavier quarks, vector and Higgs bosons as well as the t + t~ final state, nevertheless, are 
also viable options [23]. We remark that if antiprotons are produced in the hadronization 
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of final state particles [69], as in the case of quark, gluon, vector and Higgs boson channels, 
the GCE signal should be accompanied by an antiproton signal in cosmic rays. As the lat¬ 
ter is not detected, the “missing” antiproton signal appears to be in tension with the DM 
interpretation of the GCE [70, 71]. Nevertheless, this signal indeed depends strongly on the 
Galactic diffusion model, on which no consensus has yet been reached in literature [72-74], 
The missing antiproton problem can in principle be alleviated by invoking leptonic final states 
as t + t - . These, however, result in energetic electrons and positrons and therefore produce 
radio synchrotron radiation in the Galactic magnetic held. If the Galactic DM density profile 
presents a cusped core, as with the NFW profile, and the magnetic held in the interested 
region is sufficiently strong, ~10 fi G, the emitted synchrotron signal should dominate over 
the astrophysical background [70, 75, 76]. Although such a signal is not observed, our knowl¬ 
edge of the properties of the DM halo and the magnetic helds in the central region of the 
Galaxy is too limited to exclude the leptonic annihilation channel and, in the following, we 
will therefore exploit the possibility offered by the t + t~ final state. 

In the extension of the type-X 2HDM that we consider, a pair of DM particles S can 
annihilate into a pair of SM particles either via H or the SM Higgs boson h. A further 
channel, yielding four SM particle final states, is provided by the SS —> AA interaction, with 
the pseudo-scalar A decaying promptly as A —> t + t~ and/or A —> bb. As a DM particle 
lighter than the SM vector bosons provides a better fit of the GCE, 6 7 in our analysis we 
restrict the DM mass of mg accordingly and forbid heavier SM final states kinematically. 
Furthermore, below we set Asl = 0 to bar the SS AA —> SM channel, resulting in non¬ 
trivial and boosted SM final states which require a new and dedicated analysis in the context 
of GCE. We postpone the study of this channel to a future work, presenting however some 
remarks on the consequences of Asl A 0 in the last subsection. 

Under the assumptions specified above, S annihilates only via H and SM Higgs boson 
h into the dominant final states bb and t + t~ . Figure 2 shows the quantitative behaviour 
of the corresponding cross sections (av) multiplied by and averaged over the DM relative 
velocity.' The bb channel dominates for ms ~ m^/2 due to the hierarchy in the SM Yukawa 
couplings. The second peak, corresponding to m s ~ mu 1 2 is instead dominated by the t + t~ 
channel. In figure 2 we considered tp = 30, and found that the effects of this parameter 
on the annihilation cross sections are small provided that tp > 10. We have set Asl = 0, 
as explained above, and Asq = 0.004. The constraints due to the SM Higgs boson invisible 
width and direct detection experiment require that Asq < 0.006. The parameters Ai and mA 
enter the mentioned annihilation cross sections only through the decay width of h and H. 
We found that the values of the cross sections vary in a negligible way for tua E [10, 250] GeV 
and Ai E [0, 27 t/3]. 

4.2 The t + t channel 

As indicated by the grey region on the left in figure 2, fitting the GCE with the t + t~ channel 
requires a light DM particle S, ms = 9.96+Jgi GeV, and an annihilation cross section times 
velocity (av) T+T - = A x 0.337^0 04 g x 10~ 26 cm 3 s _1 , where A E [0.17,5.3] depending on 
the uncertainties in the DM halo profile at GC. 8 Within the present framework such values 
can be achieved if the S annihilation into t + t~ are mediated by H and provided that the 
mass of the latter satisfies mn ~ 2 ms- Unfortunately, we find that such a light value of 

6 See Table I and Figure 3 in ref. [23]. 

7 The expressions for the annihilation and spin independent direct detection cross sections are given in [53]. 

8 In this study we adopted the values in ref. [23] 
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m s I GeV 


Figure 2: The average annihilation cross section times the DM relative velocity for SS —> 
t + t~ (blue) and SS —> bb (red) versus the DM mass of ms- The average is performed on the 
DM velocity distribution and the mass of H is set to 23 GeV and 240 GeV in the upper and 
lower panel respectively. The remaining parameters are: tan f3 = 30, Asl = 0, AgQ = 0.004, 
rriA = 20 GeV and Ai = 0.3. The oval regions show the values of DM mass and (av) that fit 
the GCE: t + t~ channel on the left and bb channel on the right. The grey regions show the 
systematic errors due to uncertainties in the DM halo profile [23]. In the upper panel, our 
model fits both the GCE and the DM relic abundance through the t + t~ channel, satisfying 
at the same time the constraints from direct detection experiments and the SM Higgs invisible 
decay width. In the lower panel, through the bb channel, our model is able to fit the GCE as 
well as the DM relic abundance, satisfying also the constraints brought by direct detection 
experiments and by the SM Higgs invisible decay width. On top of that, in this case, our 
model lowers as well the discrepancy that afflicts the prediction of the anomalous magnetic 
moment of the muon g — 2 below the ~ 2cr. 
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run is in tension with the LEP bounds on light Higgs bosons [77, 78] and does not allow 
for a combined fit of the GCE and of the anomalous magnetic moment of the muon, as 
discussed in section 5. Nevertheless, the S annihilation into t + t~ provides a superb fit 
of the GCE as well as of the DM relic abundance. The methodology at the basis of our 
investigation is the following: we generated 400 random points within the parameter ranges 
mjj £ [15,30] GeV and Asq E [—0.005,0.005], in agreement with the constraint due to 
the invisible SM Higgs decay width and to the direct detection experiments. For every point 
(rriff- Asq) obtained in this way, we then generated other 2000 random points in the parameter 
space ms E [8.74,11.73] GeV, as required by the GCE at lcr, and tp E [1,70]. We finally 
calculated (av) and the relic. DM abundance for the resulting sample. Figure 3 shows where 
the points for which {av) T+T - and the relic DM abundance [10] are within the corresponding 
la ranges fall on the mn-^SQ and tp—ms planes. Figure 4 shows same points projected on 
the planes spanned by ms and either the spin-independent direct detection cross section crsi 
or the annihilation cross section times velocity ( av ) T+T _. 
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Figure 3: Solutions of the model for the t + t~ channel that fall within the la range (plus 
systematic uncertainties due to the DM halo profile) both of the annihilation cross section 
times velocity, required by the GCE fit, and of the DM relic abundance. The colour gradient 
denotes the value of tp E [1, 70], with darker tones corresponding to higher values. 

4.3 The bb channel 

Fitting the GCE with the bb channel yields a larger DM mass and annihilation cross section: 
ms = 48.7 l ®'2 GeV and (■ av) bb = A x 1.75+026 x 10~ 26 cm 3 s -1 with A defined above. Here 
S annihilates mainly to bb via the SM Higgs h and the required {av) bb forces ms ~ m^/2, 
as shown in the lower panel of figure 2. Adopting the bb channel to fit the GCE, the mass 
of H can be large enough to allow also for the fit of muon anomalous magnetic moment, as 
shown in section 5. On the other hand, the quality of the GCE fit is poorer than in the t + t~ 
case, remaining at about the 1.5a level. The result is obtaining by generating 800 random 
points on the parameter ranges mu £ [100,450] GeV, following the constraints of section 5 
and 3, and Asq E [—0.005,0.005]. For every point (m#, Asq) we then generated 2000 random 
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Figure 4: The same as in figure 3b, but projected on the planes spanned by the DM mass ms 
and the spin independent direct detection cross section <7si, upper panel, or the annihilation 
cross section times velocity {crv), lower panel. The coloured bands in the upper panel denote 
the constraints from the direct detection experiments LUX [27] , XENON1T [79] and LZ [80]. 
In the lower panel, the dark ellipses show the 1, 2 and 3 a confidence level around the value 
required by the GCE fit. The grey region shows the uncertainty due to the limited knowledge 
of the DM halo profile [23]. 


points covering the parameter space raj G [61, 65] GeV, within the 2cr range of the GCE best 
fit, and tp G [1, 70]. For each member of the resulting sample we calculate {crv} and the relic 
DM abundance. Figure 5 shows the points yielding a {av) hb within the corresponding 2a 
range that also reproduce the desired DM abundance [10] within la. The rectangular points 
denote the solutions for which our model’s prediction of anomalous magnetic moment of the 
muon falls within the experimental 2 a confidence level. Figure 6 shows same points projected 
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on the plane spanned ms and the spin independent direct detection cross section ergi or the 
averaged annihilation cross section times velocity (av)^. 



Asq ms / GeV 

(a) (b) 

Figure 5: Solutions of the model that fit the GCE at a 2a level (plus systematics due to DM 
halo profile uncertainties) and the DM relic abundance at lcr via the SS —> bb process. The 
colour gradient indicates the magnitude of tp £ [1,70], with larger values in a darker tone. 
The rectangular points allow for the fit of the muon anomalous magnetic moment at the 2a 
confidence level. 


4.4 Allowing for both the r + r and bb channels 

The framework we consider allows in principle for DM annihilation proceeding through both 
the t + t~ and bb channels, with a branching that can be as large as 0(1). We find however 
that having set Asl = 0 in order to avoid the SS —> A A annihilations, such a branching is 
never sizeable in our framework. The t + t~ and bb final states thus dominate the signal in 
complementary regions of the parameter space, as shown in figure 2, justifying the presented 
analysis in which the subdominant channel is disregarded. Interestingly, relaxing the condition 
Asl = 0 yields instead a sizeable branching between the considered channels and, therefore, 
a mixed t + t~ and bb final state that could be used to fit the GCE signal. Remarkably, as 
shown in figure 2, this mixed final state has the potential to improve the quality of the GCE 
fit: whereas the low energy part of the latter is well described by DM annihilation into the bb 
final state, the complementary high energy power-law tail [23, 81] is better fitted by the t + t~ 
channel. The quantitative analysis of the GCE in terms of mixed final states, as well as the 
study of the SS —> AA —» SM channel, require the development of new dedicated tools that 
goes beyond the purpose of the present paper and is therefore postponed to a future work. 

5 Impact on the anomalous magnetic moment of the muon 

We focus now on the tension between the Standard Model prediction and the measurements 
of the muon anomalous magnetic moments a^ := (g^ — 2)/2. Dedicated experiments set the 
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Figure 6: The same points as in figure 5b projected on the planes spanned by the DM 
mass ms and the spin independent direct detection cross section <tsi, upper panel, or the 
annihilation cross section (cry), lower panel. The coloured bands in the upper panel denote 
the constraints from the indicated direct detection experiments. In the lower panel, the dark 
ellipses show the 2 and 3a confidence interval for the GCE fit while the grey region indicates 
the uncertainty due to the limit knowledge of the DM halo profile [23]. The fit is always 
above the lcr level due to the condition ms ~ m^/2. 


latter to [32, 82] 

a® xp = 116592091 (54)(33) x 10“ n (5.1) 

while the corresponding Standard Model prediction, which comprises the QED, Electroweak 
and Hadronic corrections, matches [83-85] 

a® M = 116591829 (57) x 10" 11 , (5.2) 
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Figure 7: Visual fitting of the prompt 7 -ray spectrum from mixed annihilations to bb and 
t + t~ against the data of GCE. The extra annihilation to t + t~ seem to improve the quality 
of fit especially for the high energy tail of the signal. The branching ratio of the channels bb 
to t + t~ is fixed here to 0.3. 


giving rise in this way to a discrepancy between experiments and theory which is more than 
3(7 large: 


Aa M := a[f p - aJ M = 262 (85) x 10“ n . (5.3) 

Clearly, the physics beyond the Standard Model may help reducing this gap. Within the 
framework we propose, for instance, the anomalous magnetic moment of the muon receives 
additional contributions from both the extended Higgs and the dark sectors. Unfortunately, at 
the precision level achieved by current experiments, the impact of the latter is negligible given 
that the scalar held S enters the computation of Astarting at the two-loop level and that, 
differently from the case of the Barr-Zee diagrams [ 86 , 87] discussed below, the corresponding 
amplitudes do not benefit from any enhancement factor that would compensate the higher 
loop suppression. The calculation of the anomalous magnetic moment of the muon then 
proceeds here in the same fashion as in a pure lepto-specihc 2HDM, with the first contribution 
to brought at the one-loop level by the extended Higgs sector [ 88 ]: 

Sa (2HDM,U) = V (^7^(7). (5.4) 

je{h,H,A,H±} 


In the above formula rjj, := rnjj/Mj. yj^ are the rescaled Yukawa couplings of the muon, 
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presented in table 1, while the loop functions Fj are given by 


Fh,H{y) ■= f X ^ X \ dx - In (y) -\ + 0(y), (5.5) 

J 1 — x + yx- o 

o 

i 

F A (y):=[- - ,X - - ^dx ^hln(y)+ ^+0{y), (5.6) 

J 1 — x + yx~ o 

o 

^YAdx^-l + ofe). (57) 

0 

The one-loop correction brought by the extended Higgs sector is then dominated by the neutral 
scalars and pseudoscalar contributions, which interestingly have opposite signs. Notice that 
the boson h, which plays here the role of the Standard Model Higgs boson, is also included in 
the above summation. Our choice is motivated by the fact that whereas the pure Standard 
Model Higgs contribution accounted for in a® M is negligible [88], the condition in eq. (3.13), 
beside addressing the h —> AA decay width, yields small deviations of sin(/3 — a) from the 
Standard Model limit sin(/3 — a) = 1 that possibly result in a sizeable contribution to a 

On top of the one-loop corrections, the anomalous magnetic moment of the muon is also 
particularly sensitive to the Barr-Zee diagrams that, at the two-loop level, induce an effective 
coupling between the neutral scalars or the pseudoscalar and the photons. These yield [89, 90] 

H 2HDM ' BZ, = S| £ -W&W) (5-8) 

ie{h,H,A,H±} 

fe{t,b,T ,...} 


where a is the fine-structure constant, rj. := rnj / Mf. nif is the mass of a fermion /, while 
Cf, qf and y l j are respectively its colour multiplicity, electric charge and rescaled Yukawa 
coupling. The loop functions are given by 


G h ,H(y) 


G A (y) 


2x(l — x) — 1 1 (x(\ — x) 
x(l -x)-y 


In 


1 


x(l -x)-y 


In 


y 

x(l — x) 

y 


dx, 


dx. 


(5.9) 


(5.10) 


The importance of the Barr-Zee diagrams stems from the mj/m ^ enhancement of the contri¬ 
butions in eq. (5.8), which can easily overcome the extra loop suppression factor that is absent 
in the corresponding one-loop contributions. Notice also that whereas the one-loop corrections 
brought by the (pseudo)scalars of the 2HDM are (negative) positive, the corresponding two 
loop corrections have opposite signs. Then, if the contributions of the 2HDM to a^ are domi¬ 
nated by the scalar one-loop correction, it could be possible to bridge the discrepancy within 
the theoretical prediction of this quantity and the corresponding measurement. Alternatively, 
if the pseudoscalar corrections dominate and are determined by the (positive) Barr-Zee dia¬ 
gram, the model under consideration can also add on the SM prediction and effectively reduce 
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mH=m H ±=200 GeV 



rrifil GeV 

m H =m H ±=4 50 GeV 



Figure 8: The contribution to the anomalous magnetic moment of the muon brought by the 
considered model as a function of tan fi and the pseudoscalar mass tua for a fixed mu = mjj±. 
The areas shaded in dark and light grey correspond to the regions in which the discrepancy 
between theory and experiments is compatible with zero at the 68% and 95% confidence 
level respectively. The orange dashed line and the red line correspond to the 95% and 99% 
exclusion limits that lepton universality imposes [58]. 


the gap quantified in A a^. In order to investigate this possibility, we compare in figure 8 the 
latter to the total contribution brought by the model: <5a[( 2HDM ^ = <5a^ 2HDM,:L ^ + <5a|f HDM,BZ ^. 

Iu agreement with previous analyses [58, 84, 91], we find that within the considered 
framework a light pseudoscalar could in principle reconcile the theoretical prediction of a^ 
with the corresponding measurement. As shown in figure 8, if the mass of this particle 
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falls in the range [10,100] GeV and tan /3 > 35, the pseudoscalar contribution to is 

large enough to significantly reduce the discrepancy A a^. The areas shaded in dark and light 
grey respectively signal the regions where the latter is compatible with a null discrepancy at 
the 68% and 95% confidence level. The proposed picture is due to the remaining particles 
in the extended Higgs sector being much heavier than the pseudoscalar, in a way that the 
corresponding corrections to the anomalous magnetic moment are suppressed. The result 
applies on the whole range considered for the parameter mu = E [200,450] GeV. We 
also checked that the same values of the parameters result in negligible contributions to 
the anomalous magnetic moment of the electron, which remains well in agreement with its 
experimental measurement on the whole considered region of the parameter space. 

It was recently shown [58], however, that lepton universality severely constrains the 
2HDM pseudoscalar solution to the puzzle posed by the anomalous magnetic moment of 
the muon through the measurements of the leptonic r decay widths. In the lepto-specific 
2HDM, these quantities receive additional contributions from the H± decay channel, at the 
tree level, and from the remaining particles of the extended Higgs sector at higher orders. By 
computing the ratios of the r leptonic decay widths evaluated at the one-loop level, it is then 
possible to constrain the 2HDM contribution through the experimental bounds on the lepton 
universality [92], The procedure results in the 95% and 99% exclusion limit shown in figure 8 
by the orange-dashed and red line respectively, which bar a significant part of the parameter 
space associated to large tg values but weaken progressively for increasing values of and 
mn- We also remark that the current measurements of B s —> fifa disfavour pseudoscalar with 
masses mA 1$ 10 GeV [93]. The constraint is relevant if the (heavy) neutral and charged Higgs 
masses are close to the lower limits of their ranges, while it is less severe for larger values of 
these parameters. 

The 2HDM pseudoscalar contribution can then reduce, but not completely remove, the 
discrepancy between prediction and measurement of the anomalous magnetic moment of the 
muon. As shown in figure 8, lowering the latter below the 95% confidence limit requires 
a light pseudoscalar with mass 10 GeV < m.A 1$ 30 GeV and moderately large values of 
tan/3, 35 < tan/3 < 45. On top of that, alleviating the constraints that lepton universality 
casts through the leptonic r decay imposes a strongly split spectrum in the theory, with 
ttia 200 GeV < mn = m#±. Interestingly, this strict set of condition still remarkably 
allows for the fit of the GCE within the considered model. More in detail, the large value of mn 
indicated by the analyses of a^ implies that the GCE signal results, in the present framework, 
from the process SS —>• bb. Our analyses revealed that the values of tan/3 necessary to 
reproduce the detected GCE signal, figure 5b, are compatible with the requirement imposed 
by <V 

As for the possibility offered within the considered framework by a light neutral scalar H , 
by setting mA = E [100,450] GeV we find that the constraints from lepton universality 
[58] forbid a substantial contribution to Afrom a scalar as light as m h ~ 20 GeV, as 
suggested by the GCE fit via the t + t~ channel. As the corresponding exclusion exceeds 
the 3cr level, we conclude that in this case it is not possible to address the problem of the 
anomalous magnetic moment of the muon. 

6 Conclusions 

We considered the lepto-specific 2HDM augmented with a real singlet scalar S that plays 
the role of DM. While the annihilations of S into SM particles explain the GCE excess, the 
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extended Higgs sector reduces the tension between prediction and measurement of the muon 
anomalous magnetic moment. 

In this framework, the DM particles S annihilate either via the neutral scalar H or 
the SM Higgs boson h, mainly into bb and t + t~ respectively. These two possibilities cor¬ 
respond to different regions in the parameter space of the model: the bb channel requires 
ms ~ 7?7 ./j,/2 whereas the t + t~ channel dominates for ms ~ mu 1 2. In order to fit the GCE 
with the mentioned channels, we scanned over the parameter regions given by tp £ [1,70], 
m,A £ [10,250] GeV, Ai £ [0, 27T/3] and Asq £ [—0.005,0.005]. For the t + t~ channel we took 
mu £ [15,30] GeV and ms £ [8.74,11.73] GeV, while for the bb final state we considered 
mu £ [100,450] GeV and ms £ [61,65] GeV. The choice of these ranges is motivated by 
the constraints on the SM Higgs boson invisible width, by direct detection experiment and 
previous fits of GCE [23], as well as by the theoretical bounds due, for instance, to perturba¬ 
tive unitarity and vacuum stability. Beside analysing the GCE, we investigated whether it is 
possible to fit the muon anomalous magnetic moment in the considered regions of the param¬ 
eter space. In the calculation of the latter we accounted for the one-loop and the two-loop 
Barr-Zee contributions brought by the additional neutral scalar H, charged scalars and 
pseudoscalar A that appear in the considered model. 

With our method we find that: 

• The t + t~ channel fits the GCE and the DM relic abundance satisfying the constraints 
of direct detection and the Higgs invisible decay width. However, we find that because 
of the strong constraints imposed by lepton universality [58] through the measurements 
of the leptonic tau decays, it is not possible to address the puzzle of the anomalous 
magnetic moment of the muon in this setup. 

• The bb channel fits the GCE, proposing a DM candidate which respects the bounds 
imposed by cosmology, direct detection experiments and Higgs invisible decays. Com¬ 
pared to the t + t~ channel, the bb final state result in a looser fit, remaining at the 
level of ~ 1.5cr at best, mainly because of the restriction on the DM particle mass 
ms ~ mh /2 that the required annihilation cross section imposes. On the other hand, 
the setup required to fit the GCE via the bb channel remarkably allows for sizeable 
contributions to the anomalous magnetic moment of the muon that survive the strict 
constraints imposed by lepton universality. In particular, bringing the discrepancy be¬ 
tween prediction and measurements of the latter within the 95% confidence interval 
requires a light pseudoscalar with mass 10 GeV < mA 1$ 30 GeV and moderately large 
values of tan /3, 35 < tan /3 < 45. 

The next generation direct detection experiments, the LZ for example, have the potential to 
test both the proposed solutions. 

We remark that in our analyses we set Asl = 0 to bar the SS —> A A —> SM channel, 
resulting in boosted SM final states that require a dedicated analysis which we delay to a 
future work. We however argued that relaxing this condition would result in a mixed t + t~ 
and bb final state that has potential to improve the quality of the GCE fit within the presented 
framework. 

A Renormalisation group equations 

The beta functions for the gauge couplings, the top Yukawa coupling and scalar quartic 
couplings at one-loop level, defined by dgi/dp = (3 gi , where g is the renormalisation scale, are 
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given by 

16tt 2 /3 5 / = 7c/ 3 , 

16vr % = 7g 3 , 

16vr 2 /3 93 = 7gl 

IQtt 2 (3 yt = [~\g 2 - ^ g ,2 - 8 gl^j yt + ^y 3 , 

IGvr 2 /^, = - (3y 4 + 2g 2 g' 2 + 3y /4 ) — (9 g~ + 3y^)Ai + 2AX 2 + 2A§ 

+ 2 A 3 A 4 + A| + A| + 2\$ l , 

16h~(3\ 2 = - (3y 4 + 2g 2 g' 2 + 3y /4 ) — (9 g~ + 3 y / 2 )A 2 + 24A^ + 2A 2 

+ 2 A 3 A 4 + A| + A 2 + 12y 2 \2 + 2A|q, ^ 

16tt 2 P\ 3 = -(3y ! — 2g 2 g' 2 + g' 4 ) — (9y 2 + 3g /2 )Xs + 12(Ai + A 2 )A 3 + 4A§ 

+ 4(Ai + A 2 ) A 4 + 2A 2 + 2A| + 6 y 2 A 3 + 4 A 5 LA 5 Q, 

I 677 ^/3a 4 = 3g 2 g'~ — 3(3 g 2 + y / 2 )A 4 + 4(Ai + A 2 )A 4 + 8 A 3 A 4 + 4A 2 + 8 A§, 

I 671 M(3\ 5 = (—9 g 2 — 3 c/“ + 6 yf + 4Ai + 4 A 2 + 8 A 3 + 12 A 4 )As, 

167t 2 /3a s = 72Aj + 2 A| l + 2A|q, 

16vr 2 (3\ sl = —-(3y 2 + g'~)XsL + 12(Ai + 2Xs)Xsl + ( 4 A 3 + 2 A 4 )AgQ + 8 Xg L , 

16vr ~(3\ sq = —-(3y 2 + g'~)XsQ + 12(A2 + 2Ag)Asg + ( 4 A 3 + 2A4)Asq + 8A|q 

+ 6 y 2 AgQ. 

Acknowledgments 

The authors thank A. Fowlie, T. Abe, R. Sato and K. Yagyu for useful discussions. This 
work was supported by grants PUT799, PUT808, IUT23-6, CERN+, TK120, MOBILITAS 
grants MJD387, MTT 8 , MTT60 due to the European Social Fund, and by the EU through 
the ERDF CoE program. 

References 

[1] G. Jungman, M. Kamionkowski, and K. Griest, Supersymmetric dark matter , Phys.Rept. 267 
(1996) 195-373, [hep-ph/9506380], 

[2] G. Bertone, D. Hooper, and J. Silk, Particle dark matter: Evidence, candidates and constraints , 
Phys.Rept. 405 (2005) 279-390, [hep-ph/0404175]. 

[3] S. Galli, F. Iocco, G. Bertone, and A. Melchiorri, CMB constraints on Dark Matter models with 
large annihilation cross-section , Phys.Rev. D80 (2009) 023505, [arXiv: 0905. 0003]. 

[4] T. R. Slatyer, N. Padmanabhan, and D. P. Finkbeiner, CMB Constraints on WIMP 
Annihilation: Energy Absorption During the Recombination Epoch , Phys.Rev. D80 (2009) 
043526, [arXiv: 0906.1197], 

[5] G. Huetsi, A. Hektor, and M. Raidal, Constraints on leptonically annihilating Dark Matter 
from reionization and extragalactic gamma background , Astron.Astrophys. 505 (2009) 

999-1005, [arXiv: 0906. 4550], 


- 18 - 


[6] M. Cirelli, F. Iocco, and P. Panci, Constraints on Dark Matter annihilations from reionization 
and heating of the intergalactic gas, JCAP 0910 (2009) 009, [arXiv: 0907. 0719]. 

[7] G. Hutsi, J. Chluba, A. Hektor, and M. Raidal, WMAP7 and future CMB constraints on 
annihilating dark matter: implications on GeV-scale WIMPs, Astron.Astrophys. 535 (2011) 
A26, [arXiv: 1103.2766], 

[8] C. Evoli, S. Pandolfi, and A. Ferrara, CMB constraints on light dark matter candidates, 
Mon.Not.Roy.Astron.Soc. 433 (2013) 1736, [arXiv: 1210. 6845]. 

[9] M. S. Madhavacheril, N. Sehgal, and T. R. Slatyer, Current Dark Matter Annihilation 
Constraints from CMB and Low-Redshift Data, Phys.Rev. D89 (2014) 103508, 

[arXiv: 1310.3815], 

[10] Planck Collaboration, P. Ade et al., Planck 2015 results. XIII. Cosmological parameters, 
arXiv:1502.01589. 

[11] Fermi-LAT Collaboration, W. Atwood et al., The Large Area Telescope on the Fermi 
Gamma-ray Space Telescope Mission, Astrophys.J. 697 (2009) 1071-1102, [arXiv: 0902. 1089]. 

[12] L. Goodenough and D. Hooper, Possible Evidence For Dark Matter Annihilation In The Inner 
Milky Way From The Fermi Gamma Ray Space Telescope, arXiv:0910.2998. 

[13] D. Hooper and L. Goodenough, Dark Matter Annihilation in The Galactic Center As Seen by 
the Fermi Gamma Ray Space Telescope, Phys.Lett. B697 (2011) 412-428, [arXiv: 1010. 2752]. 

[14] K. N. Abazajian, The Consistency of Fermi-LAT Observations of the Galactic Center with a 
Millisecond Pulsar Population in the Central Stellar Cluster, JCAP 1103 (2011) 010, 

[arXiv: 1011.4275], 

[15] A. Boyarsky, D. Malyshev, and O. Ruchayskiy, A comment on the emission from the Galactic 
Center as seen by the Fermi telescope, Phys.Lett. B705 (2011) 165-169, [arXiv: 1012.5839]. 

[16] D. Hooper and T. Linden, On The Origin Of The Gamma Rays From The Galactic Center, 
Phys.Rev. D84 (2011) 123005, [arXiv: 1110. 0006], 

[17] K. N. Abazajian and M. Kaplinghat, Detection of a Gamma-Ray Source in the Galactic Center 
Consistent with Extended Emission from Dark Matter Annihilation and Concentrated 
Astrophysical Emission, Phys.Rev. D86 (2012) 083511, [arXiv: 1207.6047]. 

[18] C. Gordon and O. Macias, Dark Matter and Pulsar Model Constraints from Galactic Center 
Fermi-LAT Gamma Ray Observations, Phys.Rev. D88 (2013), no. 8 083521, 

[arXiv: 1306.5725], 

[19] O. Macias and C. Gordon, Contribution of cosmic rays interacting with molecular clouds to the 
Galactic Center gamma-ray excess, Phys.Rev. D89 (2014), no. 6 063515, [arXiv: 1312.6671]. 

[20] K. N. Abazajian, N. Canac, S. Horiuchi, and M. Kaplinghat, Astrophysical and Dark Matter 
Interpretations of Extended Gamma-Ray Emission from the Galactic Center, Phys.Rev. D90 
(2014), no. 2 023526, [arXiv: 1402.4090], 

[21] T. Daylan, D. P. Finkbeiner, D. Hooper, T. Linden, S. K. N. Portillo, et al., The 
Characterization of the Gamma-Ray Signal from the Central Milky Way: A Compelling Case 
for Annihilating Dark Matter, arXiv: 1402.6703. 

[22] T. Lacroix, C. Boehm, and J. Silk, Fitting the Fermi-LAT GeV excess: On the importance of 
including the propagation of electrons from dark matter, Phys.Rev. D90 (2014), no. 4 043508, 
[arXiv: 1403.1987], 

[23] F. Calore, I. Cholis, C. McCabe, and C. Weniger, A Tale of Tails: Dark Matter Interpretations 
of the Fermi GeV Excess in Light of Background Model Systematics, Phys.Rev. D91 (2015), 
no. 6 063003, [arXiv: 1411 .4647], 


19 


[24] Q. Yuan and B. Zhang, Millisecond, pulsar interpretation of the Galactic center gamma-ray 
excess, JHEAp 3-4 (2014) 1-8, [arXiv: 1404.2318], 

[25] J. Petrovic, P. D. Serpico, and G. Zaharijas, Galactic Center gamma-ray "excess" from an 
active past of the Galactic Centre?, JCAP 1410 (2014), no. 10 052, [arXiv: 1405.7928], 

[26] XENONIOO Collaboration, E. Aprile et ah, Dark Matter Results from 225 Live Days of 
XENON100 Data, Phys.Rev.Lett. 109 (2012) 181301, [arXiv: 1207.5988], 

[27] LUX Collaboration, D. Akerib et al., First results from the L UX dark matter experiment at the 
Sanford Underground Research Facility , Phys.Rev.Lett. 112 (2014) 091303, [arXiv: 1310.8214], 

[28] C. Boehm, M. J. Dolan, C. McCabe, M. Spannowsky, and C. J. Wallace, Extended gamma-ray 
emission from Coy Dark Matter, JCAP 1405 (2014) 009, [arXiv: 1401. 6458], 

[29] G. D’Ambrosio, G. Giudice, G. Isidori, and A. Strumia, Minimal flavor violation: An effective 
field theory approach, Nucl.Phys. B645 (2002) 155-187. 

[30] C. Arina, E. Del Nobile, and P. Panci, Dark Matter with Pseudo scalar-Mediated Interactions 
Explains the DAM A Signal and the Galactic Center Excess, Phys. Rev. Lett. 114 (2015) 
011301, [arXiv: 1406.5542], 

[31] A. Hektor and L. Marzola, Coy Dark Matter and the anomalous magnetic moment, Phys.Rev. 
D90 (2014), no. 5 053007, [arXiv: 1403. 3401], 

[32] Particle Data Group Collaboration, J. Beringer et al., Review of Particle Physics (RPP), 
Phys.Rev. D86 (2012) 010001. 

[33] M. Benayoun, P. David, L. DelBuono, and F. Jegerlehner, Muon g-2 estimates : Can one trust 
Effective Lagrangians and global fits ?, arXiv: 1507.02943. 

[34] J. F. Gunion, H. E. Haber, G. L. Kane, and S. Dawson, The Higgs Hunter’s Guide, Front.Phys. 
80 (2000) 1-448. 

[35] G. Branco, P. Ferreira, L. Lavoura, M. Rebelo, M. Sher, et al., Theory and phenomenology of 
two-Higgs-doublet models, Phys.Rept. 516 (2012) 1-102, [arXiv : 1106. 0034]. 

[36] CMS Collaboration, S. Chatrchyan et al., Observation of a new boson at a mass of 125 GeV 
with the CMS experiment at the LHC, Phys. Lett. B716 (2012) 30-61, [arXiv: 1207. 7235]. 

[37] ATLAS Collaboration, G. Aad et al., Observation of a new particle in the search for the 
Standard Model Higgs boson with the ATLAS detector at the LHC, Phys. Lett. B716 (2012) 
1-29, [arXiv: 1207.7214], 

[38] M. Aoki, S. Kanemura, K. Tsumura, and K. Yagyu, Models of Yukawa interaction in the two 
Higgs doublet model, and their collider phenomenology, Phys.Rev. D80 (2009) 015017, 

[arXiv: 0902. 4665], 

[39] J. F. Gunion and H. E. Haber, The CP conserving two Higgs doublet model: The Approach to 
the decoupling limit, Phys.Rev. D67 (2003) 075019, [hep-ph/0207010]. 

[40] R. Barbieri, L. J. Hall, and V. S. Rychkov, Improved naturalness with a heavy Higgs: An 
Alternative road to LHC physics, Phys.Rev. D74 (2006) 015007, [hep-ph/0603188]. 

[41] E. Ma, Verifiable radiative seesaw mechanism of neutrino mass and dark matter, Phys.Rev. 
D73 (2006) 077301, [hep-ph/0601225], 

[42] V. Silveira and A. Zee, SCALAR PHANTOMS, Phys. Lett. B161 (1985) 136. 

[43] J. McDonald, Gauge Singlet Scalars as Cold Dark Matter, Phys. Rev. D50 (1994) 3637-3649, 
[hep-ph/0702143], 

[44] C. P. Burgess, M. Pospelov, and T. ter Veldhuis, The minimal model of nonbaryonic dark 
matter: A singlet scalar, Nucl. Phys. B619 (2001) 709-728, [hep-ph/0011335]. 


- 20 - 


[45] V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf, and G. Shaughnessy, LHC 
Phenomenology of an Extended Standard Model with a Real Scalar Singlet, Phys. Rev. D77 
(2008) 035005, [arXiv: 0706. 4311], 

[46] V. Barger, P. Langacker, M. McCaskey, M. Ramsey-Musolf, and G. Shaughnessy, Complex 
Singlet Extension of the Standard Model, Phys. Rev. D79 (2009) 015018, [arXiv: 0811. 0393], 

[47] M. Gonderinger, Y. Li, H. Patel, and M. J. Ramsey-Musolf, Vacuum Stability, Perturbativity, 
and Scalar Singlet Dark Matter, JHEP 1001 (2010) 053, [arXiv:0910.3167[. 

[48] Y. Cai, X.-G. He, and B. Ren, Low Mass Dark Matter and Invisible Higgs Width In Darkon 
Models, Phys.Rev. D83 (2011) 083524, [arXiv: 1102. 1522], 

[49] C.-S. Chen and Y. Tang, Vacuum stability, neutrinos, and dark matter, JHEP 1204 (2012) 

019, [arXiv: 1202.5717], 

[50] M. Gonderinger, H. Lim, and M. J. Ramsey-Musolf, Complex Scalar Singlet Dark Matter: 
Vacuum Stability and Phenomenology, Phys.Rev. D86 (2012) 043511, [arXiv: 1202. 1316], 

[51] J. M. Cline, K. Kainulainen, P. Scott, and C. Weniger, Update on scalar singlet dark matter, 
Phys.Rev. D88 (2013) 055025, [arXiv: 1306.4710], 

[52] M. Boucenna and S. Profumo, Direct and Indirect Singlet Scalar Dark Matter Detection in the 
Lepton-Specific two-Higgs-doublet Model, Phys.Rev. D84 (2011) 055011, [arXiv: 1106. 3368], 

[53] N. Okada and O. Seto, Gamma ray emission in Fermi bubbles and Higgs portal dark matter, 
Phys.Rev. D89 (2014), no. 4 043525, [arXiv: 1310.5991], 

[54] C. Bonilla, D. Sokolowska, J. L. Diaz-Cruz, M. Krawczyk, and N. Darvishi, IDMS: Inert Dark 
Matter Model with a complex singlet, arXiv: 1412.8730. 

[55] A. Alves, S. Profumo, F. S. Queiroz, and W. Shepherd, Effective field theory approach to the 
Galactic Center gamma-ray excess, Phys. Rev. D90 (2014), no. 11 115003, [arXiv: 1403.5027], 

[56] G. Belanger, K. Kannike, A. Pukhov, and M. Raidal, Minimal semi-annihilating Zjv scalar 
dark matter, JCAP 1406 (2014) 021, [arXiv: 1403.4960], 

[57] K. Kannike, G. Belanger, A. Pukhov, and M. Raidal, Semi-annihilating scalar dark matter, 

PoS EPS-HEP2013 (2013) 400. 

[58] T. Abe, R. Sato, and K. Yagyu, Lepton-Specific Two Higgs Doublet Model as a Solution of 
Muon g — 2 Anomaly, arXiv: 1504.07059. 

[59] ATLAS, CMS Collaboration, G. Aad et al., Combined Measurement of the Higgs Boson Mass 
in pp Collisions at s = 7 and 8 TeV with the ATLAS and CMS Experiments, Phys. Rev. Lett. 
114 (2015) 191803, [arXiv: 1503. 07589], 

[60] S. Kanemura, T. Kubota, and E. Takasugi, Lee-Quigg-Thacker bounds for Higgs boson masses 
in a two doublet model, Phys.Lett. B313 (1993) 155-160, [hep-ph/9303263]. 

[61] A. G. Akeroyd, A. Arhrib, and E.-M. Naimi, Note on tree level unitarity in the general two 
Higgs doublet model, Phys.Lett. B490 (2000) 119-124, [hep-ph/0006035]. 

[62] I. Ginzburg and I. Ivanov, Tree level unitarity constraints in the 2HDM with CP violation, 
hep-ph/0312374. 

[63] I. Ginzburg and I. Ivanov, Tree-level unitarity constraints in the most general 2HDM, 

Phys.Rev. D72 (2005) 115010, [hep-ph/0508020]. 

[64] K. Kannike, Vacuum Stability Conditions From Copositivity Criteria, Eur.Phys.J. C72 (2012) 
2093, [arXiv: 1205.3781], 

[65] F. Lyonnet, I. Schienbein, F. Staub, and A. Wingerter, PyR@TE: Renormalization Group 
Equations for General Gauge Theories, Comput.Phys.Commun. 185 (2014) 1130-1152, 

[arXiv: 1309.7030], 


- 21 


[66] M. E. Peskin and T. Takeuchi, Estimation of oblique electroweak corrections , Phys.Rev. D46 
(1992) 381-409. 

[67] Gfitter Group Collaboration, M. Baak et al., The global electroweak fit at NNLO and 
prospects for the LHC and ILC, Eur.Phys.J. C74 (2014) 3046. [arXiv: 1407.3792]. 

[68] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, et al., Exotic decays of the 125 GeV Higgs 
boson , Phys.Rev. D90 (2014), no. 7 075004, [arXiv: 1312.4992]. 

[69] M. Cirelli, G. Corcella, A. Hektor, G. Hutsi, M. Kadastik, et al., PPPC 4 DM ID: A Poor 
Particle Physicist Cookbook for Dark Matter Indirect Detection, JCAP 1103 (2011) 051, 

[arXiv: 1012.4515], 

[70] T. Bringmann, M. Vollmann, and C. Weniger, Updated cosmic-ray and radio constraints on 
light dark matter: Implications for the GeV gamma-ray excess at the Galactic center, Phys.Rev. 
D90 (2014), no. 12 123001, [arXiv: 1406.6027], 

[71] M. Cirelli, D. Gaggero, G. Giesen, M. Taoso, and A. Urbano, Antiproton constraints on the 
GeV gamma-ray excess: a comprehensive analysis , JCAP 1412 (2014), no. 12 045, 

[arXiv: 1407.2173], 

[72] C. Evoli, I. Cholis, D. Grasso, L. Maccione, and P. Ullio, Antiprotons from dark matter 
annihilation in the Galaxy: astrophysical uncertainties, Phys.Rev. D85 (2012) 123511, 

[arXiv: 1108.0664], 

[73] G. Giesen, M. Boudaud, Y. Genolini, V. Poulin, M. Cirelli, et al., AMS-02 antiprotons, at last! 
Secondary astrophysical component and immediate implications for Dark Matter, 

arXiv:1504.04276. 

[74] C. Evoli, D. Gaggero, and D. Grasso, Secondary antiprotons as a Galactic Dark Matter probe, 
arXiv:1504.05175. 

[75] G. Bertone, M. Cirelli, A. Strumia, and M. Taoso, Gamma-ray and radio tests of the e+e- 
excess from DM annihilations, JCAP 0903 (2009) 009, [arXiv:0811 .3744]. 

[76] R. Crocker, N. Bell, C. Balazs, and D. Jones, Radio and gamma-ray constraints on dark matter 
annihilation in the Galactic center, Phys.Rev. D81 (2010) 063516, [arXiv: 1002.0229]. 

[77] J. Cao, P. Wan, L. Wu, and J. M. Yang, Lepton-Specific Two-Higgs Doublet Model: 
Experimental Constraints and Implication on Higgs Phenomenology, Phys.Rev. D80 (2009) 
071701, [arXiv: 0909. 5148], 

[78] S. Chang, S. K. Kang, J.-P. Lee, and J. Song, Higgs potential and hidden light Higgs scenario 
in two Higgs doublet models, arXiv: 1507.03618. 

[79] E. Aprile and XENON1T collaboration, The XENON1T Dark Matter Search Experiment, 
ArXiv e-prints (June, 2012) [arXiv: 1206. 6288]. 

[80] D. C. Mailing, Akerib, et al., After LUX: The LZ Program, ArXiv e-prints (Oct., 2011) 

[arXiv: 1110.0103], 

[81] F. Calore, I. Cholis, and C. Weniger, Background model systematics for the Fermi GeV excess, 
JCAP 1503 (2015) 038, [arXiv: 1409.0042], 

[82] Muon g-2 Collaboration, G. Bennett et al., Final Report of the Muon E821 Anomalous 
Magnetic Moment Measurement at BNL, Phys.Rev. D73 (2006) 072003, [hep-ex/0602035]. 

[83] K. Hagiwara, R. Liao, A. D. Martin, D. Nomura, and T. Teubner, {g — 2) mu and a(M f) 
re-evaluated using new precise data, J.Phys. G38 (2011) 085003, [arXiv: 1105. 3149]. 

[84] A. Broggio, E. J. Chun, M. Passera, K. M. Patel, and S. K. Vempati, Limiting 
two-Higgs-doublet models, JHEP 1411 (2014) 058, [arXiv: 1409.3199]. 

[85] M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, Reevaluation of the Hadronic Contributions 
to the Muon g-2 and to alpha(MZ), Eur.Phys.J. C71 (2011) 1515, [arXiv: 1010.4180]. 


- 22 - 


[86] Y.-L. Wu and Y.-F. Zhou, Muon anomalous magnetic moment in the standard model with two 
Higgs doublets , Phys.Rev. D64 (2001) 115018, [hep-ph/0104056]. 

[87] S. M. Barr and A. Zee, Electric dipole moment of the electron and of the neutron, Phys. Rev. 
Lett. 65 (Jul, 1990) 21-24. 

[88] F. Jegerlehner and A. Nyffeler, The Muon g-2 , Phys.Rept.. 477 (2009) 1-110, 

[arXiv :0902. 3360] . 

[89] D. Chang, W.-F. Chang, C.-H. Chou, and W.-Y. Keung, Large two loop contributions to g-2 
from a generic pseudoscalar boson, Phys.Rev. D63 (2001) 091301, [hep-ph/0009292], 

[90] K. Cheung and O. C. Kong, Can the two Higgs doublet model survive the constraint from the 
muon anomalous magnetic moment as suggested?, Phys.Rev. D68 (2003) 053003, 
[hep-ph/03021 1 1] . 

[91] L. Wang and X.-F. Han, A light pseudoscalar of 2HDM confronted with muon g-2 and 
experimental constraints, JHEP 1505 (2015) 039, [arXiv: 1412.4874]. 

[92] Heavy Flavor Averaging Group (HFAG) Collaboration, Y. Amhis et ah, Averages of 
b-hadron, c-hadron, and t- lepton properties as of summer 2014, arXiv: 1412.7515. 

[93] H. E. Logan and U. Nierste, B(s, d) —> l + (,~ in a two Higgs doublet model, Nucl.Phys. B586 
(2000) 39-55, [hep-ph/0004139]. 


- 23 


